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ABSTRACT 



We present a model atmosphere analysis of cool hydrogen-rich white dwarfs observed at the Very Large Telescope (VLT) with the 
X-shooter spectrograph. The intermediate-dispersion and high signal-to-noise ratio of the spectra allowed us to conduct a detailed 
analysis of hydrogen and heavy element line profiles. In particular, we tested various prescriptions for hydrogen Balmer line broad- 
ening parameters and determined the effective temperature and surface gravity of each star. Three objects (NLTT 1675, 6390 and 
1 1393) show the presence of heavy elements (Mg, Al, Ca, or Fe). Our abundance analysis revealed a relatively high iron to calcium 
ratio in NLTT 1675 and NLTT 6390. We also present an analysis of spectropolarimetric data obtained at the VLT using the focal 
reducer and low dispersion spectrograph (FORS) and we established strict upper limits on the magnetic field strengths in three of 
the DAZ white dwarfs and determined the longitudinal magnetic field strength in the DAZ NLTT 10480. The class of DAZ white 
dwarfs comprises objects that are possibly accreting material from their immediate circumstellar environment and the present study 
contributes in establishing class properties. 
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1. Introduction 

Studies of heavy elements in cool hydrogen-rich (DA) white 
dwarfs were originally r estricted to a f ew c ases, i.e., the 
cool white dwarf G74-7 dLacombe et al.L I1983I) and the ZZ 
Ceti star G29-38 dKoester et all Il997l) . because of their rar- 
ity. Diffusion timescales in cool DA white dwarfs are rela- 
tively short dPaquette et all I1986I) . and heavy elements were 
not expected to be detectable in the photosphere of these ob- 
jects. Therefore, heavy elements such as calcium in G74-7 must 
be continuously suppli e d from an external source. Incidentally, 
IZuckerman & Becklinl dl987 | ) rep orted an infrared excess in 
G29-38, and iGraham et alJ d 1990ft showed that the excess is 
caused by the presence of circumstellar dust. The detection of 
heavy elements (Mg , Ca, Fe) in the atmosphere of G29-38 by 
iKoester et alJ d 19971) strengthened the link between atmospheric 
contaminants and the presence of circumstellar material. 

Surveys of white dwarfs conducted using high-dispersion 
spectrographs attached to large aperture teles copes (e.g., 
IZuckerman et all 1 1 998L 120031: IKoester et all 120051) reveal that, 
contrary to earlier assessments, a significant fraction of DA 
white dwarfs, labelled DAZ white dwarfs, are contaminated 
with elements heavier than helium. Within a sample of ~ 
100 DA white dwarfs that are not in close binary systems, 
Zuckerman et al. (2003) showed that w 25% of the stars exhibit 
photospheric heavy element lines; for the majority of stars in 
their sample only calcium lines are observed, but in some cases 



* Based on observations collected at the European Organisation for 
Astronomical Research in the Southern Hemisphere, Chile under pro- 
gramme ID 080.D-0521, 082.D-0750, 083.D-0540, 084.D-0862 and 
086.D-0562. 



other elements such as magnesium, iron, aluminium and silicon 
are also detected. 

Our study of cool white d warfs in the revised NLTT cat- 
alogue dSalim & Gould 120031) unveiled new cases of external 
contaminations. The stars were selected using a reduced proper 
motion diagram co mbined with optical and infrared colours 
dKawka et all 12004) resulting in a catalogue of ~ 400 objects; 
only about half of these objects had previously been spectro- 
scopically confirmed as white dwarfs. Although the NLTT cat- 
alogue is deemed incomplete at low Galactic latitude and far 
south (see iLepine & Sharal [2005), it still contains many stars 
that remain largely unstudied. More complete catalogues utiliz- 
ing the Digital Sky Surveys, such as the LSPM catalo gue of stars 
with proper motions greater than 0.15" yr~' (Lepine & Shara, 
2005) have been compiled and complement the NLTT catalogue 
where it is incomplete. 

iKawka & Vennesl ((2006) conducted low-dispersion spectro- 
scopic observations of white dwarf candidates from the NLTT 
catalogue and listed 49 new objects including 30 DA white 
dwarfs, three DAZ white dwarfs including NLTT 43806, and 
16 non-DA white dwa rfs. In following-up on NLTT 43806, 
Zuckerma n et al.l d201 lh obtained high signal-to-noise and high- 
dispersion spectra showing it to harbour a weak magnetic field. 

In this work, we present intermediate-dispersion and high 
signal-to-noise (S/N) ratio spectroscopy and a model atmosphere 
analysis of new DAZ white dwarfs that were selected from the 
New Luyten Two-Tenths (NLTT) catalogue. These new objects 
add to the few hydrogen-rich white dwarfs that have been scru- 
tinized for heavy elements other than calcium. We presented a 
first report on this programme w ith the analysis of the co olest 
DAZ in our sample (NLTT 10480: Kawka & Vennesl l201lh . The 



1 



A. Kawka and S. Vennes: Chemical composition of cool white dwarfs 
Table 1. Log of spectroscopic observations 6 



NLTT 


Instrument 


UT Date 


Exp. Time 

(s) 


1675 


NTT/EFOSC2 


2009 Aug. 24 


2 x 1500 




VLT/X-shooter 


2010 Nov. 13 


2 x 2400 






2010 Dec. 8 


2 x 2400 


6390 


VLT/FORS1 


2007 Nov. 1 


2 x 1500 




VLT/X-shooter 


2011 Jan. 1 


1 x 2400 


11393 


VLT/FORS1 


2007 Nov. 3 


2 x 1260 




NTT/EFOSC2 


2008 Oct. 21 


2 x 1200 




VLT/X-shooter 


2010 Nov. 13 


2 x 2400 






2011 Jan. 11 


2 x 2400 


23966 


CTIO/R-C Spec 


2008 Feb. 25 


2 x 2400 




VLT/FORS2 


2010 Jan. 23 


4x 1200 




VLT/X-shooter 


2011 Jan. 6 


1 x 1800 






2011 Jan. 30 


1 x 1800 






2011 Mar. 26 


1 x 1800 



X-shooter spectrum of NLTT 10480 showed Zeeman-split cal- 
cium and hydrogen lines that revealed a surface magnetic field 
of B s = 0.519 + 0.004 MG. At the present time, onl y three other 
magn etic DAZ white dwarfs are known: G77- 50 dFarihi etldl 
1201 ll). NLTT 43806 dZuckerman et all 120111) and LTT 8381 
dKoester et all l2009h . Apart from these four magnetic DAZ 
white dwarfs, two more white dwarfs with helium dominated 
atmospheres that are contaminated with heavy elements have 
also b een shown to b e weakly magneti c: G165-7 dDufour et all 
2006) and LHS 2534 dReid et allBOOTl) . Recently. IPotter & Toutl 
(2010) and lNordhaus" et alJ (1201 ll) have proposed that magnetic 
fields in white dwarfs may be generated during the common- 
envelope phase of i nteracting binaries. Applying this model to 
the case of G77-50. [Farihi et aD d201 ll) proposed that the pres- 
ence of a weak field in this star may be the result of a past 
common-envelope episode with a planetary component during 
the white dwarf's formative years. 

In the following Sect. 2 we describe observations obtained 
at the European Southern Observatories (ESO) using the New 
Technology Telescope (3.6-m) and the Very Large Telescopes 
(VLTs) as well as the 4-m telescope at Cerro Tololo Inter- 
American Observatory (CTIO). Sect. 3 presents our model at- 
mosphere analysis including details of the model structure (Sect. 
3.1) and heavy-element line profile (Sect. 3.2) calculations. 
Using these models we determine the effective temperature and 
surface gravity for each star (Sect. 3.3), and measure their heavy 
element abundance (Sect. 3.4). The stellar radial velocities and 
kinematics, and estimates of the magnetic field strengths are pre- 
sented in Sect. 3.5 and Sect. 3.6, respectively. We summarize and 
discuss some implications of our results in Sect. 4. 

2. Observations 

2.1. Spectroscopy 

We first observed objects from the present sample during a pro- 
gramme aimed at identifying new white dwarfs in the NLTT cat- 
alogue. NLTT 6390 and NLTT 1 1393 were observed with the fo- 
cal reducer and low dispersion spectrograph (FORS1) attached 
to the 8m UT2 (Kuenyen) at Paranal Observatory on UT 2007 
November 1 and 3 as part of our spectropolarimetric survey of 
white dwarfs. The purpose of these observations was to confirm 
the nature of the observed sample of stars and to search for weak 
magnetic fields. We used the 600B grism combined with a slit- 
width of 1 arcsecond that provided a resolution of 6 A and a 
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A(A) 

Fig. 1. Identification spectra, fx versus A, normalized at A = 
5500A and shifted up for clarity. 



spectral range of 3780 to 6180 A. The observations consisted of 
a sequence of two consecutive exposures. In the first exposure 
the Wollaston prism is rotated to -45° followed by a second ex- 
posure with the Wollaston prism rotated to +45° from which we 
extracted the flux and circular polarization spectra. 

As part of our identification programme we observed 
NLTT 23966 with the R-C spectrograph attached to the 4m tele- 
scope at CTIO on UT 2008 February 25. We used the KPGL2 
(316 lines per mm) grating with the WG360 filter to block out 
the second order contamination. We set the slit- width to 1.5 arc- 
seconds which provided a spectral resolution of ~ 8 A and a 
spectral range of 3820 to 7500 A. The spectra were flux cali- 
brated with the flux standard GD 108. We also selected a sam- 
ple of spectroscopically confirmed hydrogen-rich white dwarfs 
from the NLTT catalogue, including NLTT 23966, for spec- 
tropolarimetric observations with the aim of searching for weak 
magnetic fields. NLTT 23966 was observed with the FORS2 at- 
tached to the 8m UT1 (Antu) at Paranal Observatory on UT 2010 
Jan 23. We used the 1200R grism with the order blocking filter 
GG435 and 2kx4k MIT CCD. The slit-width was set to 1 arc- 
second which resulted in a spectral resolution of ~ 3 A and a 
spectral range of 5810 to 7290 A. The observations were car- 
ried out in the same way as for NLTT 6390 and NLTT 11393 
with the FORS 1 spectrograph with one sequence consisting of 
two consecutive exposures with the Wollaston prism rotated 90° 
between the two exposures. Two sequences were obtained for 
NLTT 23966. 

NLTT 1675 was first observed with the ESO Faint Object 
Spectrograph and Camera (EFOSC2) attached to the 3.6m New 
Technology Telescope (NTT) at La Silla Observatory on UT 
2009 August 24. We used Grism 11 which has 300 lines per 
mm and a blaze wavelength of 4000 A. The slit-width was set 
to 1 arcsecond which resulted in a spectral resolution of ~ 14 A 
and a spectral range of 3700 to 7250 A. We have also observed 
NLTT 1 1393 using this setup on UT 2008 October 21. The ob- 
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Table 2. Line identifications 



1. \ 



1 1 I 




I 1 1 1 I 
NLTT 23966 



NLTT 6390 



NLTT 1675 





NLTT 11393 




3720 3740 3760 3940 3960 3980 

A(A) 

Fig. 2. X-shooter spectra covering iron {left) and calcium (right) 
lines in four cool, hydrogen-rich white dwarfs from the NLTT 
catalogue. Photospheric calcium lines are observed in three stars, 
while iron is only observed in NLTT 1 675 and 6390. Contrary to 
our initial report (Kawk aet al.L 120171) . our new X-shooter spec- 
tra show that NLTT 23966 is devoid of calcium. The possible 
identification of aluminium in reported for NLTT 1675. 



servations were carried out at the parallactic angle and were flux 
calibrated with the flux standard Feige 110. 

Figure Q] shows the low-dispersion identification spectra of 
the five cool hydrogen rich white dwarfs for which higher disper- 
sion spectra were obtained with VLT/X-shooter, including the 
identification spectra o f the cool magneti c DAZ NLTT 10480 
which was analyzed by Kawka & Vennesl d201 lb . The low dis- 
persion spectra of these stars showed the presence of the 
Caii/13933A line and therefore were selected for follow-up 
high-dispersion observations with X-shooter. In the case of 
NLTT 23966 the Canvl3933A line seen in the CTIO spectrum 
was the result of background noise. 

Following-up on our low-dispersion spectroscopic obser- 
vations, we obtained a set of echelle spectra for all five ob- 
jects (including NLTT 10480) using the X-shooter spectro- 
raph attached to t he UT2 (Kueyen) at Paranal Observatory 
Vernet et"aQ 1201 ll) . The spectra were obtained with the slit- 
width set to 0.5, 0.9 and 0.6 arcsecond for the UVB, VIS and 
NIR arms, respectively. This set-up delivered a resolving power 
of 9100 for UVB, 8800 for VIS and 6200 for NIR. 

Figure [2] shows X-shooter spectra of NLTT 1675, 6390, 
11393 and 23966. Photospheric lines of calcium, iron and alu- 
minium are marked. Table [TJ summarizes our spectroscopic ob- 
servations and Table|2]lists the line identifications and equivalent 
width measurements in the X-shooter spectra. The observations 
of NL TT 10480 were originally presented in iKawka & Vennesl 
(1201 ll) . 



Ion 


A 

(A) 




E.W. 
(mA) 










NLTT 








1675 


6390 


11393 


Fei 


3440.606 


148 


141 




Fei 


3565.379 


38 : 


48 




Fei 


3570.097 


80 


109 




Fei 


3581.195 


185 


161 




Fei 


3608.859 


84 : 






Fei 


3618.768 


97 : 






Fei 


3631.463 


104 : 


95 




Fei 


3647.842 




46 




Fei 


3719.935 


117 


191 




Fei 


3734.864 


150 


119 




Fei 


3737.131 


128 


105 




Fei 


3745.561 


112 


108 




Fei 


3748.262 




50 




Fei 


3749.485 


71 : 


106 




Fei 


3758.233 


53 


72 




Fei 


3798.511 


46 






Fei 


3815.840 


42 : 


29 : 




Fei 


3820.425 


115 


117 




Fei 


3825.880 


46 


63 




Mgi 


3832.304 


96 : 


60 : 




Fei 


3834.222 


42 






Mgi 


3838.292 


58 : 


91 : 




Fei 


3859.911 


119 


77 




Can 


3933.663 


1161 


485 


420 


All 


j 944. 006 


48 : 






Ah 


3961.520 


36 : 






Can 1 


3968.469 


866 


499 


420 


Fei 


4045.813 


56 






Cai 


4226.728 


73 






Fei 


4383.544 


51 


31 




Can 


8542.09 


64 






Can 


8662.14 


45 






(1) Blended with weak He. 



Table 3. Photometric and astrometric properties of the sample 



NLTT 


V 1 


J 1 


cos S 3 


Us 


3 




(mag) 


(mag) 


(mas yr~') 


(mas yr~ 


) 


1675 4 


17.73+0.04 


16.56 ±0.13 


203 ±6 


-100 ± 


6 


6390 


17.33 ±0.04 


16.36 ±0.13 


440 ±6 


-213 ± 


6 


11393 


17.16 ±0.04 


16.04 ± 0.08 


250 ±6 


-239 ± 


6 


23966 


17.13 ±0.04 


16.54 ±0.11 


-160 ±6 


-51 ± 


6 



(1) This work; V-magnitude for NLTT 1675 estimated from SDSS 
ugriz photometry (2) 2M ASS, ISkrutskie et alJd2006l) . <3) rNLTT, 
ISalim & Gould i2q03h. (4) Also known as LSPM J003 1 +221 8, 
iLepine & Sharal ( [20051) . 



2.2. Photometry 

We searched for photometric measurements using VizieR. The 
Two Micron All Sky Survey (2MASS) listed infrared JHK mag- 
nitudes, but only the J magnitude was of acceptable quality. For 
NLTT 6390, 11393 and 23966 we used the acquisition images 
from the X-shooter observations to estimate a V magnitude. The 
X-shooter acquisition images for NLTT 1675 were unusable. To 
set the zero point for these magnitudes we used 1 1 acquisition 
images of Feige 110 obtained between UT 2010 Dec 10 and 
2011 Jan 1. We employed the atmospheric extinction table of 
iPatatetal] (120111) . 
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For NLTT 1675, we also attempted to measure a V magni- 
tude from the EFOSC2 acquisition images, however both acqui- 
sition images appear to have been taken through thin clouds and 
the final magnitudes are unreliable. Fortunately, Sloan Digital 
Sky Survey (SDSS) photometry is available for NLTT 1675 and 
we have used these measurements to determine a V magnitude: 
We obtained a best fit model spectrum to SDSS photo metry and 
then we convolved this model with the V bandpass dBesselll 
1990) providing us with an estimate of the V magnitude. 

Table [3] lists the photometric and proper motion measure- 
ments for each star. 

3. Analysis 

3.1. Model atmospheres and spectral syntheses 

We extended our grid of model atmospheres for cool hyd rogen- 
rich white dwarfs employed in iKawka & Vennesl (I201 ll) to en- 
compass the effective temperature range 4900 < T e g < 8 000 K 
in 100K steps and the surface gravity range 7.0 < logg < 8.75 
in steps of 0.25 dex. The models are in convective and radiative 
equilibrium. Again, all relevant species (H, H + , H2, H£ , H3 ) are 
included in the statistical equilibrium equation, with electrons 
contributed by identifiable trace elements (e.g., calcium) in- 
cluded in the charge conservation equation. However, in warmer 
models that are relevant to the present study (T e g > 6 000) 
electrons are contributed mostly by the ionization of hydro- 
gen atoms. The H2-H and H-H collision-induced absorptions in 
the far Lya wing (see [kowalski & Saumonl | 2006l) a re included 
using opacity tables from iRohrmann et al.l (I201 lh . Synthetic 
colours as well as detailed hydrogen and heavy element line 
profiles are computed using the model structures. Table lA.il lists 
some photometric properties of these models. The colour indices 
at shorter wavelengths are effected by the Lya collision-induced 
absorptions. 

The calculation of synthetic hydrogen line profiles deserves 
further attent i on. The FW HM for resonant levels computed by 
lAli & Grieml (fl%5l fl966h . the AG formalism, assuming "reso- 
nance" interaction (oc R^) is 

r 3B = -= 4.8624 xl0 8 (^V /2 ^, (1) 

where, for exa mple, the lower (I) and upp er (u) levels were taken 
as Is and 2p in lAli & Grierr] (ll965L[T966l) . with gi/g u = 1/3 and 
the oscillator strength //„ = 0.4162, providing a FWHM for the 
2p level of 2.4 x 10~ 8 rad s'' cm 3 . 

iBergeron et al.l (Il99lh proposed to apply this formalism to 
principal quantum numbers / = 1 and u by taking gi/g u = (l/u) 2 
and //„ as the total oscillator strength between principal quantum 
numbers I = 1 and u. This formalism offers a means to estimate 
the total width of the level u that includes the s, p, d contri- 
butions, i.e., resonant and non-resonant levels. Figure [3] shows 
Balmer line widths calculated by summing the lower (u) and up- 
per (m') level widths F3 = r 3j „ + r 3i „/ . In this formalism, the total 
width is dominated by the u = 2 width, while the contribution of 
the upper level u' to the width decreases with increasing values 
of u'. 

A comparison wit h the calculat ions of Barkle m et al.l d2000a) 
for Ha, y6, and y, and lAllard et al.l d2008l) for Ha shows that the 
total line width is not well represented by this generalisation of 
the AG formalism. Figure [3] shows that the two fo rmalisms be- 
have d i fferen tly along the line series. Although the lAli & Grieml 
(fl%5l fl966h values are nearly constant (when expre ssed in 
rads -1 ), the width computed by Bark lem et al.l d2000ah triples 
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Fig. 3. Self-broadening paramete r T for members of t he Ba lmer 
line series. The widths (T 3 ) from lAli & Grieml d 19651 1 19661) fol- 
low a m arkedly different trend than widths from Barklem et al. 
(2000a). Improvements are obtained by adding a van der Waals 
term (r 6 ) to lAli & Grieml (fT965l[T966l) : T 3 + T 6 . 



between Ha and y. The lAli & Gri em (1965, 1966) calculations 
are dominated by 2p level width and are based o n the resonance 
ls-2p t ransition, while the widths tabulated by iBarklem et al.l 
( 2000a]) rely on the np-md formalism for the n — m transition 
and their calculations employ the "dispersive-inductive" (van der 
Waals) interaction (oc R~ 6 ) in addition to the "resonance" inter- 
action. 

Consequently, the AG formalism may be improved by adopt- 
ing a van de Waals term of the form (see lKurucz & Avreti.ll98l|) 

F 6 «8xl0- 9 r 4 3/I0 (« 4 -/ 4 ) 2/5 , (2) 

where T4 = 77 10 4 . T he total line width is the n obtained by 
adding the resonance (lAli & Grieml 1 19651 Il966h and van der 
Waals terms: r = r 3 + F(,. Figure [3] shows that the resulting line 
wi dth increases a l ong the series although the widths calculated 
by Barkle m et al.l d2000al) remain considerably larger. 

IBarklem et al.l (|2000a) also determined validity limits for the 
quasi-static framework: 35, 13, and 8 A for Ha, /3 and y, respec- 
tively, at T — 4665 K, revealing potential difficulties with higher 
members of the Balmer line series. 

In order to assess possible systematic errors in the mea- 
surement of stellar parameters based on Balmer line profiles, 
we explored three different approaches for the calculation of 
Balmer line spectra. In the first approach, the Balmer line pro- 
files from Ha to the serie s limit are computed using widths from 
lAli& Grieml (119651 119661) alone, Fag, while in the second ap- 
proach we combined AG and van der Waals widths, Fag.vcTw- In 
the last approach, t he line profiles from H a to Hy are computed 
using widths from Bar klem et al.l (l2000al) . Tepo, multiplied by 
factors of 0.75 and 1.0, again to assess uncertainties, while the 
upper Balmer lines are computed using TacvcIw- The Voigt pro- 
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files including Doppler and se lf-broadening are convolved with 
Stark-broadened profiles from lLemkel {T997). 



3.2. Heavy element line profiles 

The dominant broadening mechanism at ~ 6000 K is collision 
with hydrogen atoms, but the effect of Stark broadening is also 
included as it contributes significantly at ~ 8000 K. 

For broadening of Mg i, Al i, Ca r, and Ca n lines due to col- 
lisions with h y drogen atoms we employed the coefficients of 
Barklemetal. (2000b) where the full- width at half-maximum 
(FWHM) of the Lorentzian profiles is given by: 



w 



n(U I) 



\10 4 K/ 



(l-aO/2 



F rads 



1 cm 3 . 



(3) 



For example, we adopted logF = -7.562 for Ca 1/14226 
and -7.76 for CaH&K at T = 10 000 K, and 
a = 0.238 for Ca 1,14226 and 0.223 for CaH&K. Stark 
broadening wi dths for Cai and Can li ne s wer e ob- 
tained from | bimitriievic & SahaU Brechot d 19991) and 
iDimitriievic & Sahal-Brechol dl992l) . respectively. Additional 
parameters (Stark and c ollisions with hydrog en) were obtained 
from the compilation of lKurucz & Belli (1995). 

Hydrogen molecules provide $20% of the gas pressure in 
some layers of models at r e ff =6000 K and a negligible fraction 
at 8000 K. Consequently, the effect of collision with molecules 
is not considered further. 
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Fig. 4. Balmer line fits of X-shooter spectra using the line- 
broadening prescription TacvcIW (see text). With the exception 
of NLTT 23966, the He/Ca H blend is excluded from the line 
fitting. 



3.3. Effective temperature and surface gravity measurements 

Fitting the V - J colour index to synthetic colours (Table IA. 1 b 
constrains the effective temperature of the cool white dwarfs 
(Table @}. Within error bars the cooler objects are at J e ft ~ 
6000 K, while NLTT 23966 is hotter with r eff * 7600 K. These 
preliminary colour-based estimates are corroborated by our de- 
tailed Balmer line analysis. 

Next, we simultaneously constrained T e ff and logg by fit- 
ting the X-shooter spectra with the model grids using x 2 mini- 
mization techniques. All lines up to the series limit, H8/9 for the 
cooler stars and H10 for NLTT 23966, are included in the analy- 
sis, but we excluded the He/Ca H blend in the three cooler stars. 
Figure [4] shows the Balmer line profiles of the four white dwarfs 
compared to the best-fitting model spectra that were calculated 
using F A G,vdw- 

Figure [5] illustrates the range of effective temperatures and 
surface gravities attained for the sample when employing dif- 
ferent broadening parameters. Values connected with thin grey 
lines were obtained using either Fag alone (upper points), 
or TacvcIw (lower points). The values connected with full 
black lines were obtained using Fbpo (Ha, J3,y) combined with 
TACvdw for the upper Balmer lines (lower points), or reduced 
0.75 x Tbpo values combined with TAG.vdw (upper points). The 
derived parameters are compa red to evolutionary cooling tracks 
dBenvenut o & Althaus, 1999) allowing estimates of the mass 
and cooling age of each object. 

The average mass for the three cooler stars computed us- 
ing FAG.vdw line widths is 0.53 M e . Similarly, the average mass 
is 0.57 M Q using the combination Fbpo X 0.75/rAG,vdw- On the 
other hand, the average mass of 0.96 M obtained using Fag ap- 
pears improbably high, while the average mass of 0.44 M ob- 
tained using the combination FBPo/TAG.vdw appears improbably 
low. Higher (lower) values of the surface gravity (hence density) 
are derived as a compensation for using smaller (larger) values 



Table 4. Stellar properties 



NLTT 


T, S (V 


-J) 




logg 1 


Mass 




(K) 




(K) 


(cm s" 2 ) 


(M s ) 


1675 


5800 ± 


300 


6020 + 50 


8.04 ± 0.07 


0.61 ± 0.04 


6390 


6200 ± 


300 


6040 ± 40 


7.90 ± 0.07 


0.53 ± 0.04 


11393 


5800 ± 


300 


5890 ± 30 


7.86 ± 0.06 


0.51 ±0.03 


23966 


7600 ± 


300 


7790 ± 30 


8.20 ± 0.04 


0.72 ± 0.03 



Average of the parameters determined using Balmer line profiles 
computed using r AG-v dw an d r BP0 x 0.75 models. 



of r. Considering the uncertainties in F, it seems approp r iate to 
adopt an average mass near 0.6 M (see iTremblav et al., 2010) 
for the three cooler stars. 

On the other hand, the surface gravity, hence mass, of 
NLTT 23966 is only marginally affected by variations in the 
self-broadening F value because of the dominant effect of Stark 
broadening at higher temperatures. The higher than average 
mass estimated for NLTT 23966 may also be attributed to the 
appro ximate treatment of co nvection in current model calcula- 
tions dTremblav et al.Ll201 lh . 

Table [4] lists the adopted effective temperature and surface 
gravity of the cool white dwarfs. As a compromise, we adopted 
the average of the values obtained using model spectra calcu- 
lated with TAG.vdw or Tbpq x0.75. Individual white dwa rf masses 
were calculated following lBenvenuto & Althaus! (t 1999). 

3.4. Abundance of heavy elements 

The abundances were computed at logg = 8.0 and T e ff 
6000 + 100 K for the cooler stars, and at logg = 8.25 and 
r eff = 7800 + 100 K for NLTT 23966. Table lists the mea- 
surements obtained by fitting the spectra with models of varying 
composition; the quoted errors include the (mild) effect of vary- 
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Fig. 5. Range of parameters (T e g vs. log g) determined for a sam- 
ple of four cool hydrogen-rich white dwarfs using three differ- 
ent prescriptions for hydrogen line broadening (see text). The 
measured parameters a re com pared to evolutionary tracks from 
iBenvenuto & Althausl JT999) with masses ranging from 0.3 to 
1.0 M , and ages ranging from 1 to 5 Gyr. An age discontinuity 
is apparent ~ 0.45 M because of the different core composition 
adopted for models with M < 0.4 M (He) and M > 0.5 M 
(C/O). 

Table 5. Abundance measurements 1 



NLTT 


logMg/H 


logAl/H 


logCa/H 


logFe/H 


1675 


-8.56±0.12 


-9.28±0.17 


-9.53±0.03 


-8.63+0.13 


6390 


-8.66±0.20 


< -9.2 


-10.00+0.04 


-8.57+0.11 


11393 


< -8.7 


< -9.4 


-10.24+0.04 


< -9.4 


23966 






< -10.4 





(I) Error bars are l<r and the upper limits are given at 99% confidence. 



ing the temperature by +100 K, and the upper limits are taken at 
99% certainty (2.6<x). However, the magnesium and aluminium 
abundances are based on ~ 2<x detections and higher signal-to- 
noise ratio spectra are required to help improve the measure- 
ments. On the other hand, a strict upper limit on the calcium 
abundance in NLTT 23966 shows that its atmosphere is devoid 
of heavy elements. Figure [6] shows the various metal lines de- 
tected in NLTT 1675 compared to the best fitting model spectra. 



Fig. 6. X-shooter spectra of NLTT 1675 compared to the best fit 
model. The wavelength scale is centred on the values indicated 
on the right. The top two spectra show the weak Al i lines, fol- 
lowed by spectra showing Mg i lines around 3838 A and 3832 
A and Fe i lines around 3740 A and 3720 A. Finally, the bottom 
three spectra show the Ca i and Ca n lines. 

Table 6. Age, distance and kinematics 



NLTT 


Age 


My 


d 


U, V, w 






(Gyr) 


(mag) 


(pc) 


(kms~') 




1675 


2.2-2.9 


14.27+0.11 


49±3 


-15±6,-44±8, 


-4±6 


6390 


1.7-2.3 


14.06+0.11 


45±2 


-37+5,-101+9, 


-25±8 


11393 


1.8-2.2 


14.13±0.10 


40±2 


2±5,-72+7, - 


5±7 


23966 


1.6-1.9 


13.45+0.06 


54±2 


-23+5,-31+6, 


-6±6 



Table [6] lists the cooling age, absolute magnitude, distance 
and Galactic space velocities of the sampl e. The ages were com- 
puted using the evolutionary models of IBenvenuto & Althausl 
(Il999h and the absolute magni tudes were calculate d usin g 
our model grid and radii from IBenvenuto & A lthaus (1999). 
The photometric distances were calculated using the apparent 
(TableO and absolute (Tabled Johnson V magnitudes. Finally, 
the Galactic velocity component s UVW were calculated us- 
ing iJohnson & Soderbloml d 1987b . The lag i n V suggest that 
the white dwarfs belon g to the old thin disk (ISion et all | 1988; 
iPauliet al.Ll200ll2006h . 



3.5. Radial velocity measurements, kinematics, and ages 

We measured radial velocities v r = 15.0,83.2,64.2 and 
66.4km s" 1 for NLTT 1675, 6390, 11393 and 23966, respec- 
tively. We estimated the error to be 5.0km s _1 based on the scat- 
ter in individual line measurements and the expected precision of 
a R ~ 9000 spectrograph. Taking into account the expected grav- 
itational redshift of the white dwarfs (based on the spectroscop- 
ically determined parameters), the actual velocities of the stars 
are v = -16.5 + 5.8, 58.3 + 5.8,40.9 + 5.4 and 25.4 + 5.4km s" 1 
for NLTT 1675, 6390, 1 1393 and 23966, respectively. 



3.6. Magnetic field strengths 

A lack of splitting at the instrument resolution (R ~ 9000) im- 
poses a limit on the magnetic field strength. Setting a limit of 
approximately 1/3 of a resolution element to the putative mag- 
netic broadening of narrow line cores due to Zeeman splitting: 

M = kB^~ (4) 

where k - 4.67 x 10~ 13 /l 2 g e ff, ^ is the wavelength in A, B s is the 
average surface field in G, and g e ff is the effective Lande factor 
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Table 7. Magnetic field measurements 
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Fig. 7. Circular polarization (top) and flux (bottom) spectra of 
NLTT 10480. The observed circular polarization spectrum is 
compared to a model polarization spectrum at Bi — 212 kG (top). 



(= t and | for Ca K and H, respectively). For Ca K the resulting 
limit on the field strength is Z? S5 17 kG. However, this effect 
could be confused with the effect of stellar rotation. 

We obtained VLT/FORS spectropolarimetry for three stars 
(NLTT 6390, 1 1393 and 23966) in addition to NLTT 10480. We 
calculated the mean longitudi nal magnetic (Bi in G) field using 
the weak- field approximation (Ang el et all 1 19731) 



kB, 



vF 



dF/dA 



(5) 



where v = V/I is the degree of circular polarization, F = I 
is the total spectral flux, and dF/dA is the flux gradient. First, 
we fitted the observed Balmer line profiles with model spectra 
to determine the best fit line profiles. We then used these line 
profiles to calculate the flux gradient and determine the longitu- 
dinal magnetic field. To measure the longitudinal field strength 
for NLTT 6390 and NLTT 11393 we used HB to H6 and for 
NLTT 23966 we used Ha. 

Table [7] lists the longitudinal field measurements for three 
stars (NLTT 6390, 11393 and 23966 ) from this paper and fo r 
the cool white dwarf NLTT 10480 dKawka & Vennesl 1201 lb . 
The quoted errors are lcr. Using the Balmer lines, no fields 
stronger than ~ 5 to ~ 20 kG are detected in NLTT 6390, 
1 1393 and 23966, but we measured a longitudinal field strength 
B, = -200.4 + 124.0 kG in NLTT 10480, which is barely a 2cr 
detection. The large uncertainty in Bi is mostly due to the weak- 
ness of the Balmer lines (HB, y and 8) in NLTT 10480, but a 
measurement based on the calcium lines appears more promis- 
ing. 

Equation includes only the first term of a Taylor expan- 
sion that is used to calculate v. This is valid where the Zeeman 
splitting is small compared to the line width. For Ca H&K in 
NLTT 10480, the Zeeman splitting is fully resolved and there- 



NLTT 


Bi 


Lines 




(kG) 




6390 


-4.8 + 9.0 


HB, Hy, H<5 


10480 


-200 ± 124 


HB, Hy 




-212 + 50 


CaH&K 


11393 


+7.8 ± 12.2 


HB, Hy, H6 


23966 


-2.6 + 3.1 


Ha 



fore we need t o include the second term of the Taylor expansion 
as described inlMathvs & Stenflol d 1986b 



1 / dF 1 , d 3 F\ 



(6) 



B y fitting the CaH&K c ircular polarization spectrum depicted 
in IKawka & Vennesl d201 lb we measured a longitudinal field of 
-212 + 50 kG in NLTT 10480. Figure|7]shows the best fit model 
circular polarization and flux spectra to the observed FORS 
spectra of NLTT 1 0480. A reasonable match is achieved between 
the model and observed spectra, except for the left cr component 
of the Call /13968 A line, which appears stronger than the model. 
This is an improvement over the Balmer line measurements, but 
it still has a significantly larger uncertainty than the measure- 
ments obtained for the other white dwarfs in the sample. This is 
mainly due to the lower signal-to-noise ratio in the spectral re- 
gion containing the Can lines. Combining the measu rement of 
the surface field of 5 19 ± 4 kG (IKawka & Vennesl 1201 lb and the 
Ca H&K measurement of the longitudinal field in NLTT 10480 
suggests an inclination below 59°, in agreement with an angle of 
60 + 3° obtained fr om the strength of the ca lcium Zeeman-split 
lines determined in lKawka & Vennesl (1201 lb . 

4. Summary and discussion 

We have presented a detailed model atmosphere analysis of 
high-quality spectroscopy of a sample of cool DA white dwarfs. 
We show that the atmospheres of NLTT 1675, NLTT 6390, 
and NLTT 11393 are contaminated with heavy elements. On 
the other hand we show that, contrary to the DAZ white dwarf 
NLTT 10480, none of these three new DAZ white dwarfs har- 
bour a magnetic field. We also refined our analysis of the mag- 
netic field in NLTT 10 480 and confirm our original results 
(IKawka & Vennesl 1201 lb . 

Figure [8] shows the abundances of iron and calcium relative 
to hydrogen for a sample of cool DAZ white dwarfs, including 
those from the present sample, and the corresponding abundance 
ratio of iron with respect to calcium (Fe/Ca ). The published 
abund ance measurements were taken fro m IZuckerman et al.l 
(12001 and Ikawka et all d201lb. and from iFarihi et al I d201 lb. 
IKawka & Venne s (2011) and Zucker man et al.l d201 lb for the 
three magnetic white dwarfs, G77-50, NLTT 10480 and 
NLTT 43806, respectively. Table IB. fl lists the selected measure- 
ments. The Fe/Ca abundance ratio appears to be constant as 
a function of temperature with an average of <Fe/Ca>= 13.1 
and a dispersion of cr = 0.33 dex. The linear fit to the data 
exclud es NLTT 43806, which has a Fe/Ca abundance ratio of 
* 1.3 (IZuckerman et al.ll201 lb making it significantly iron poor 
compared to the rest of the sample. Similarly, Fe/Al is » 0.6 in 
NLTT 43806 while it is * 4.5 in NLTT 1675. Both calcium and 
iron abundances show an increase as a function of temperature, 
but the scatter of iron abundance measurements (cr ~ 0.3 1 dex) 
is lower than that of calcium (cr ~ 0.95 dex). It is possible that 
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Fig. 8. Abundance ratio of iron with respect to calcium versus 
the effective temperature (top). The best fitting line that excludes 
NLTT 43806 is also shown. The abundance of iron and cal- 
cium versus the effective temperature (bottom), upper limits are 
not shown. The best fitting lines for iron (full line) and calcium 
(dashed line) are shown. 

Table 8. Inferred accretion rates log Mx(g s -1 ) 



NLTT 


Mg 


Al 


Ca 


Fe 


1675 


6.59 


5.97 


5.98 


7.17 


6390 


6.46 


<6.02 


5.48 


7.20 


11393 


<6.36 


<5.76 


5.18 


<6.31 


23966 






<5.08 





higher signal-to-noise and higher dispersion spectra of objects 
apparently devoid of iron would uncover even lower iron abun- 
dances and thereby increase the scatter in the measurements. 
In summary, the Fe/Ca abundance fraction is ~10 compared to 
ss 14 in the Sun. The Mg/Ca abundance fraction of NLTT 1675, 
NLTT 6390, G77-50 and NLTT 43806 range between * 6 and 
* 30. 

We calculated the presumed accretion rates assuming steady 
state between accretion and diffusion for magnesium, alu- 
minium, calcium and iron. We used the diffusion timescales 
and mass of the convective zone of lKoesterl <|2009). For alu- 
minium we used the d iffusion timescale log tai = 4.23 that 
IZuckerman et alJ (120111) adopted for NLTT 43806, which is a 
star with temperature and surface gravity similar to our sam- 
ple of DAZ white dwarfs. Table [8] lists the accretion rates for 
the NLTT 1675, 6390 and 11393. For objects with upper limits 
on abundances, such as NLTT 23966, we list the corresponding 
upper limits on accretion rates. The observed abundance ratios 
relative to calcium in meteorites (CI chondrites) or in the Sun 
dGrevesse et al.L 120071) are log Mg/Ca, log Al/Ca, log Fe/Ca = 
1.23, 0.10, and 1.15, respectively. Our sample appears to con- 
form with these ratios. Following sample trends, magnesium and 
iron may well be present in NLTT 1 1393 but slightly below de- 
tection limits in our spectra. 



Considering the whole sample, the observed calcium and 
iron abundances follow the excepted temperature trend as- 
suming a constant acc retion rate over cooling ages (see 
iKoester & Wilkenl l2006). although it appears somewhat steeper. 
The calcium and iron abundances increase by 0.32 and 0.4 dex 
per 1000 K interval, respectively, compared to predicted in- 
creases of 0.26 and 0.27 dex per 1000 K interval assuming uni- 
form accretion over the sample. This steeper trend suggests that 
the accretion flow onto white dwarfs may be slightly abating 
with time, or that the actual particle flux at the bottom of the con- 
vection zone may be underestimated in cooler, i.e., older white 
dwarfs. 

Prior to this study, five DAZ white dwarfs with T e ff ;S 
6000 K) were known where two of the se, NLTT 43806 
dKawka & Vennesl 120061) and NLTT 10480 dKawka & VennesL 
1201 lh were discovered as part of our spectroscopic observa- 
tions of white dwarf candidates from the rNLTT catalogue. 
In this paper we add three more cool DAZ white dwarfs to 
this important sample of accreting objects. Also notable is the 
high incidence of magnetism among this small sample of stars: 
Three out of the eight DAZ white dwarfs cooler than « 6000 K 
(G77-50, NLTT 43806 and NLTT 10480) are weakly magnetic. 
This incidence is markedly higher th an previously estimated 
dKawka et all 120071: lLiebertltaLll2003l) . 
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Table A.l. Selected synthetic colours 



logs T eB (K) B-V V-J g-J g-r r-i 



(K) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


5 800 


0.535 


1.149 


1.428 


0.396 


0.141 


6000 


0.496 


1.070 


1.325 


0.354 


0.119 


6200 


0.463 


0.998 


1.232 


0.316 


0.100 


6400 


0.432 


0.931 


1.145 


0.281 


0.081 


6 600 


0.405 


0.867 


1.063 


0.248 


0.063 


6 800 


0.380 


0.807 


0.986 


0.218 


0.047 


7 000 


0.358 


0.750 


0.913 


0.190 


0.031 


7 200 


0.337 


0.695 


0.845 


0.164 


0.017 


7 400 


0.319 


0.642 


0.779 


0.140 


0.003 


7 600 


0.303 


0.591 


0.716 


0.117 


-0.011 


7 800 


0.289 


0.542 


0.656 


0.095 


-0.024 


8 000 


0.276 


0.494 


0.598 


0.075 


-0.037 


5 800 


0.537 


1.139 


1.419 


0.395 


0.139 


6000 


0.496 


1.061 


1.316 


0.352 


0.118 


6 200 


0.463 


0.992 


1.226 


0.314 


0.100 


6400 


0.432 


0.926 


1.140 


0.280 


0.082 


6 600 


0.405 


0.863 


1.060 


0.247 


0.065 


6 800 


0.380 


0.804 


0.984 


0.218 


0.049 


7 000 


0.357 


0.748 


0.912 


0.190 


0.034 


7 200 


0.337 


0.694 


0.845 


0.164 


0.019 


7 400 


0.319 


0.643 


0.781 


0.141 


0.006 


7 600 


0.303 


0.593 


0.720 


0.118 


-0.007 


7 800 


0.289 


0.545 


0.662 


0.098 


-0.020 


8 000 


0.276 


0.499 


0.606 


0.078 


-0.032 


5 800 


0.542 


1.128 


1.409 


0.396 


0.136 


6000 


0.497 


1.050 


1.304 


0.350 


0.116 


6 200 


0.461 


0.982 


1.215 


0.311 


0.098 


6400 


0.430 


0.918 


1.131 


0.277 


0.081 


6600 


0.403 


0.857 


1.053 


0.245 


0.065 


6 800 


0.378 


0.799 


0.978 


0.215 


0.050 


7 000 


0.355 


0.744 


0.908 


0.188 


0.035 


7 200 


0.335 


0.691 


0.842 


0.162 


0.021 


7 400 


0.316 


0.641 


0.779 


0.139 


0.008 


7 600 


0.300 


0.593 


0.720 


0.117 


-0.004 


7 800 


0.286 


0.546 


0.663 


0.097 


-0.016 


8 000 


0.274 


0.501 


0.609 


0.078 


-0.028 
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A. Kawka and S. Vennes: Chemical composition of cool white dwarfs 
Table B.l. Abundance measurements of cool white dwarfs. 



WD 


Name 


T eS (K) 


log(Ca/H) 


log (Fe/H) 


Reference 


0028+220 


NLTT1675 


6010 


-9.53 


-8.63 


1 


0032-175 


G266-135 


9235 


-10.20 




2 


0151-308 


NLTT6390 


6020 


-10.00 


-8.57 


1 


0208+396 


G74-7 


7201 


-8.84 


-7.89 


2 


0243-026 


LHS1442 


6798 


-9.90 




2 


0245+541 


G174-14 


5190 


-11.70 




2.3 


0315-293° 


NLTT 10480 


5200 


-10.30 




4 


0322-019° 


G77-50 


5310 


-9.8 


-8.70 


5 


0334-224 


NLTT 11393 


5880 


-10.24 




1 


0543+579 


GD290 


8142 


-10.29 




2 


0846+346 


GD96 


7373 


-9.41 




2 


1054-226 


NLTT25792 


8160 


-7.30 




6 


1102-183 


EC11023-1821 


8026 


-10.43 




2 


1124-293 


ESO439-80 


9680 


-8.53 




2 


1202-232 


EC 12028-23 16 


8619 


-9.78 


-8.00 


2.6 


1208+576 


G 197-47 


5830 


-10.96 




2 


1225+006 


HE1225+0038 


9390 


-9.8 




6 


1257+278 


G149-28 


8600 


-8.04 


-7.41 


7 


1315-110 


HE1315-1105 


9080 


-9.4 




6 


1344+106 


G63-54 


6945 


-11.13 




2 


1407+425 


PG 


9856 


-9.87 




2 


1 /ICC . Tflf) 

1455+298 


LHS3007 


7366 


-9.31 


-8.60 


2.6 


1633+433 


G 180-63 


6569 


-8.63 


-7.55 


2 


1653+385° 


NLTT43806 


5900 


-7.9 


-7.80 


7 


1821-131 


LHS3384 


7029 


-10.71 




2 


1826-045 


G21-16 


9480 


-8.83 




2 


1858+393 


G205-52 


9470 


-7.84 


-6.92 


2 


2115-560 


LTT8452 


9570 


-8.0 


-6.90 


6 


2221-165 


HE2221-1630 


9990 


-8.0 
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Appendix B: Abundance measurements 



n 



